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Abstract. Let (M,g,k) be an initial data set for the Einstein equations of general relativity. 

We prove that there exist solutions of the Plateau problem for marginally outer trapped surfaces (MOTSs) 
that are stable in the sense of MOTSs. This answers a question of G. Galloway and N. O'Murchadha raised 
in |12| and is an ingredient in the proof of the spacetime positive mass theorem j9] given by L.-H. Huang, 
D. Lee, R. Schoen and the first author. 

We show that a canonical solution of the Jang equation exists in the complement of the union of all 
weakly future outer trapped regions in the initial data set with respect to a given end, provided that this 
complement contains no weakly past outer trapped regions. The graph of this solution relates the area of 
the horizon to the global geometry of the initial data set in a non-trivial way. 

We prove the existence of a Scherk-type solution of the Jang equation outside the union of all weakly 
future or past outer trapped regions in the initial data set. This result is a natural exterior analogue for the 
Jang equation of the classical Jenkins-Serrin theory. 

We extend and complement existence theorems 1 161 1171 35 25 15 27 for Scherk— type constant mean 
curvature graphs over polygonal domains in (M, g), where (M, g) is a complete Riemannian surface. We can 
dispense with the a priori assumptions that a sub solution exists and that [M, g) has particular symmetries. 
Also, our method generalizes to higher dimensions. 



1. Introduction and Notation 

In this paper we prove a series of results related to the geometric theory of the Jang equation and appli- 
cations to classical problems in the theory of minimal and constant mean curvature graphs. 

In Section [2] we show that there exist stable solutions of the Plateau problem for marginally outer trapped 
surfaces, answering a question of G. Galloway and N. O'Murchadha raised in [12]. This result is subtle in 
view of the non-variational nature of these surfaces. The proof is based on the fact that using the existence 
theory in [7] we can construct ordered families of solutions of the Plateau problem. This result is used in the 
recent proof of the spacetime positive mass theorem by L.-H. Huang, D. Lee, R. Schoen and the first author 
in 0. 

In Sections [3HU and [5] we develop the geometric theory of the Jang equation pioneered by R. Schoen and 
S.-T. Yau in |30j to prove the existence of non-trivial and, in some cases, canonical Scherk-type solutions of 
the Jang equation in the complements of the total weakly future outer trapped region and the total weakly 
past outer trapped region. 

In Section [6] we employ techniques from the geometric theory of the Jang equation, in particular the 
capillarity regularization and the geometric blow up analysis from |30j and ideas from the solution of the 
non-variational Plateau problem for marginally outer trapped surfaces in [J , to the classical Jenkins-Serrin 
problem [161 117] of finding necessary and sufficient conditions for a domain in a Riemannian surface to sup- 
port a Scherk-type constant mean curvature graph. In the case of positive mean curvature, we are able to 
dispense with the a priori assumption that the domain admit a sub solution which is required in the foun- 
dational paper by J. Spruck [35j (in K 2 ) and its recent extension to domains in § 2 and H 2 by L. Hauswirth, 
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H. Rosenberg, and J. Spruck |15j . Moreover, our results are valid in arbitrary complete Riemannian sur- 
faces. In particular, the existence of a Scherk-type graph is curiously reduced to a (generically) finite set 
of inequalities relating area and circumference of certain polygons that can be inscribed into the domain: 
the Jenkins-Serrin-Spruck flux conditions. Our approach here does not distinguish between minimal and 
(positive) constant mean curvature graphs. In the case of minimal graphs, we recover the results by B. Nelli 
and H. Rosenberg [25 (in H 2 ) and by A. Pinheiro [27] (in arbitrary Riemannian surfaces). Our methods 
carry over to higher dimensions. A detailed overview of the literature and a precise statement of our result 
are given in Section [6] 

In the appendices we collect several results that are needed in the paper. The simple indirect proof of 
the (known) interior gradient estimate for solutions of the prescribed mean curvature equation based on 
the regularity theory for almost minimal boundaries in low dimensions in Appendix does not seem to 
appear in the literature, surprisingly. In Appendix [B] we characterize the boundary of domains that support 
infinite boundary value solutions of the prescribed mean curvature equation without making an a priori 
assumption regarding the regularity of the boundary. In Appendix [C] we state another simple and useful 
consequence of the classical compactness and regularity theory for almost minimal boundaries: the horizontal 
parts of the unit normal vector fields of non-parametric solutions to the prescribed mean curvature equations 
are equicontinuous in low dimensions. This is an important ingredient for the geometric analysis in Section[B] 

We proceed by introducing the notation employed throughout this paper, and by reminding the reader of 
some background material. 

Let (Af, g) be a connected Riemannian manifold of dimension n, with 3 < n < 7, and let k be a symmetric 
(0, 2)-tensor on M. In the context of the Cauchy problem for the Einstein equations in general relativity, k is 
referred to as the (spacetime) second fundamental form tensor and the triple (M, g, k) is called an initial data 
set. We often require that (M,g, k) is asymptotically flat, i.e. that the complement of some compact subset 
of M consists of finitely many connected components N\, . . . , N m , called the ends, each one diffeomorphic 
to M n \Bi(0) and such that in the corresponding coordinate systems the metric tensor g^ converges to the 
Euclidean metric 5ij and the second fundamental form tensor kij to zero. More precisely, we require that 

n — 2 

\gij - Sij\ + \x\\d k gij\ = 0(\x\~ q ) and % = 0(|x|~ 9_1 ) as |af| oo, for some q > . 

When n — 3, we ask in addition that for some /3 > 2, 

ti g (k) = g lj hj = Od^l - ' 3 ) as x oo. 

This last condition is imposed so that certain barriers for the Jang equation can be constructed far out in 
the asymptotically flat ends, cf. Appendix [Dl 

Let (M, g,k) be an initial data set, and let E C Af be a two-sided hypersurface with unit-normal vector 
field v. The future outer and past outer expansion scalars of E are defined, respectively, as 9^ = H^+trs(k) 
and = — trs(fc). Here, is the mean curvature scalar of E computed as the tangential divergence of 
v and trs(fc) is the trace of k restricted to the tangent space of E. The hypersurface E is called future outer 
trapped (respectively past outer trapped) if 9^. < (9^ < 0) everywhere on E, weakly future outer trapped 
(respectively weakly past outer trapped) if 9^ < (9^ < 0) everywhere on E, and is called a marginally 
future outer trapped surface or MOTS for short (respectively marginally past outer trapped surface or MITS 
for short) if 0+ = (9^ — 0) everywhere on E. Except for Section [5J the MOTSs and MITSs appearing in 
this paper will be closed and in fact boundaries of sets f2 that contain part of an end {x € iVj : \x\ > ro} 
for some ro > I and i G {I, . . . , m}. If we write a hypersurface E as the (relative) boundary of a set f2, say 
E = dfl in U where U is an open subset of M, we will use the unit normal field v of E pointing into f2 to 
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compute the scalar mean curvature Hy, and the expansion scalars 9^ unless otherwise noted. 

Let (M, g, k) be a complete asymptotically flat manifold of dimension 3 < n < 7. Fix one of the ends, 
say N\. It is easy to see that Hs r > \ tr$ r k\ for all r > ro, provided that ro > 1 is sufficiently large. Here, 
S r := {x £ Nx ■ \x\ — r} is the coordinate sphere of radius r in Nx and the mean curvature scalar Hs r is 
computed as the tangential divergence of the unit normal pointing into the end. Let M_ C M (respectively 
M+ C M) be the interior of the intersection of all open subsets SlcM that contain {x £ Nx : |x| > ro} 
and which have smooth properly embedded boundary satisfying Hgn + tran(k) < (Hgu — tran(fe) < 0). 
It follows from [3] in dimension n = 3 and from [8] in dimensions 3 < n < 7 that M_ (respectively M+) 
has compact embedded boundary and that Hom_ + tr^M_ (k) = (Hqm + — t*dM + (k) = 0). Thus dM- is a 
MOTS and dM+ is a MITS. The complements of the regions M_ and M+ are the total weakly future outer 
trapped region and the total weakly past outer trapped region of (M, g, k) with respect to the chosen end, 
respectively. If (M,g, k) has more than one end, then both 9M_ and dM + are non-empty and each of them 
separates the portion {x € Nx '■ \x\ > r } of Nx from \J™ =2 { X e : \ x \ — r o}- 

Given an open subset Q C M and a function u £ C 2 (Q) we define, in local coordinates near a point x £ Q, 

H{u) := U _ - and tr(Jb)(») = - 

Here, all geometric operations (raising indices, gradient, length of gradient, Hessian) are with respect to g. 
These definitions are independent of the particular coordinate system used. The function H{u) is the scalar 
mean curvature of the graph G = {(x, u(x)) : x £ ft} of u in the Riemannian product (fi x R, g + (dx n+1 ) 2 ) at 
the point (x, u{x)) with respect to the downward pointing unit normal, and tr(fc) is the trace of k (extended 
to S] x R by zero in the vertical direction) over the tangent space of this graph at (x, u(x)). If 

(1) Jf(u)+tr(fc)(u)=0, 

then G, with its downward orientation, is a MOTS in the new initial data set (M 
Equation ([T]) is known as the Jang equation. 

For background material on MOTSs, MITSs, and the Jang equation we refer the reader to the survey 
article pQ. 
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2. Stability of solutions of the Plateau problem 

The following definition of stability for MOTSs spanning a non-empty boundary is the natural analogue 
of the notion of stability for closed MOTSs that has been introduced and studied systematically in [2]. 

Definition 2.1 (Cf. [121 Section 2]). Let (M,g, k) be an initial data set and consider a two-sided hypersurface 
ScM with boundary <9£ = T and with a designated "outward" unit normal vector field v . Assume that S 
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is a MOTS. Then E is said to be stable in the sense of MOTSs if there exists a smooth function f that is 
positive in the interior o/E, vanishes on the boundary T, and is such that L^f > 0. Here, 

Lx<t> := -A s + 2(X, D^4>) + (± R s -\\h + fc|| - J(u) - /x + div s X - \X\ 2 ) 0, 

where X is the tangential part of the vector field dual to k(y,-) on E, h is the second fundamental form 
of Yi (with trace H^j, is the tangential gradient of <j) along E, Rs is the scalar curvature of E, /i = 

■j (Km +(^M{k)) 2 — I^Im) is the local mass density, and where J := div(fc — trM{k)g) is the local current 
density. Equivalently, E is stable in the sense of MOTSs if, and only if, the principal eigenvalue of L-^ is 
non-negative. 

For a careful discussion of principal eigenvalues of (not necessarily self-adjoint) elliptic operators and their 
eigenfunctions, we refer the reader to [351 Sections 3.6 and 3.7]. 

We recall the following theorem from [7J: 

Theorem 2.2 ([7 j). Let (M, g, k) be a complete initial data set of dimension n with 2 < n < 7. Let VI C M 
be a bounded open set with smooth boundary dVl. Let T C dVl be a non-empty smooth closed embedded 
submanifold of DVl such that dVl\ T = d-VlUd+Vl for disjoint non-empty relatively open subsets <9_f2, <9+f2 of 
dVl. Assume that H^n + trgo k < near d-Vl with the mean curvature computed as the tangential divergence 
pointing into VI and that Uga +tign k > near d+VL with the mean curvature scalar computed as the tan- 
gential divergence of the unit normal pointing out of VI. Then there exists a smooth hyper surf ace SCSI with 
boundary T that is an almost minimizing relative boundary in Q and such that E is a MOTS with respect to 
the unit normal pointing towards d+Vl. 

The following theorem answers a question posed in [121 Section 3]. It is an ingredient in the proof of the 
spacetime positive mass theorem given in [5]. 

Theorem 2.3. Assumptions as in Theorem \2.'2[ Then there exists a solution E of the Plateau problem for 
MOTSs in VI with boundary T that is stable in the sense of MOTSs. 

Proof. Given e > small, let T £ := {6 e : distan(0, T) = e}. From the explicit construction in the 

proof of Theorem 12.21 in pj Chapter 4] we see that the MOTSs E e C VI spanning V s are (strictly) ordered. 
To see this, recall that the open subset of Vt whose relative boundary is E E is the geometric limit as t \ 
of downward translations of the regions lying above the graphs u\ £ C^ c (Vt) constructed in [3 Lemma 4.2]. 
Given t > and < e < e' small, we have that 5^' C S^i (in the notation of [7]) and hence uf < ttf. It 

follows that the regions above the graphs are ordered so that E e lies to one side (towards d+il) of E £ . The 
geometric maximum principle shows that components of E e and E e that span components of T E and T e 
cannot touch unless they coincide. We will discard all extraneous closed components of E e . This does not 
change that each E e is a relative boundary, nor that the E e 's are ordered. 

The geometric compactness properties of the almost minimizing relative boundaries E e show that as 
e \ 0, the E £ converge smoothly and with multiplicity one to an embedded MOTS E that spans T. We 
claim that this MOTS E is stable in the sense of MOTSs. To see this, let U, V, W C E be non-empty open 
subsets with smooth boundaries such that U <e V <s W <s int E. Let v be the unit normal vector field 
of E that points towards d+Vt. (This makes sense because E is a relative boundary in VI spanning T.) By 
assumption, there exist positive functions f £ £ C°°(W) for e > small with {exp e (f e v) : 9 £ V} C E e and 
such that f £ — > with all derivatives on compact subsets of W. Because E and E 6 both satisfy the MOTS 
equation, f £ is solution of a homogenous linear elliptic equation in divergence form on V. The operator 
describing the linearization of the equation at the function that vanishes identically is is- Arguing exactly 
as in 34, p. 333], using Harnack theory, it follows that the functions f E can be rescaled (so their infimum 
is one on V, say) so as to converge smoothly to a positive function fy £ C°°(V) with Ls(fv) — 0. This 
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implies that U is stable in the sense of MOTSs. To see this, let A be the first principal eigenvalue of L-z\u 
and let h £ C°°(U) be the corresponding first (Dirichlet) eigenfunction so that L-gh — Xh. We recall that 
the first principal eigenvalue is simple and that the corresponding eigenfunctions do not change signs. By 
scaling, using that g vanishes on the boundary of U and that fy is positive on V and that U <s V, we may 
assume that < h < fy on U with equality at some point. The maximum principle then implies that A > 0. 
We conclude that every open subset U (e int £ is stable in the sense of MOTSs. Using that the principal 
eigenvalue of an elliptic operator depends continuously on the operator and the domain, it follows that £ is 
stable in the sense of MOTSs. □ 

3. SCHERK-TYPE SOLUTIONS OF THE JANG EQUATION 

The content of the following proposition is similar to that of Theorem 3.1 in [24]. We include an alternative 
proof here as preparation for the more general and difficult Theorem 13.31 below. The modification of the 
data near the boundary in our proof is much less delicate than that in [24] . The regions M_ and M+ in the 
statement of Proposition 13.11 and Theorem 13.31 below are defined in the introduction. 

Proposition 3.1 (Cf. Theorem 3.1 in [23] when n = 3). Let (M,g,k) be a complete asymptotically flat 
initial data set of dimension n, 3 < n < 7, and assume that <9Af_ and dM + are disjoint. There exists a 
smooth solution u : M_ n M+ — > R of the Jang equation H(u) + tr(fc) (it) = such that u(x) — > as x — > oo 
in the end of (M,g, k) contained in M— l~l M+ and such that u(x) — > ±oo as dist(x, dM±) — > 0. 

Proof. We abbreviate Q := M_ n M+. Let \ € Cf (M) be such that \ = ±1 near 9M±. Given e £ (0, 1) we 
define k £ := k + e\- Note that Hqm + + trgM + (k e ) > and Hqm_ + tr^Af_ (k e ) < 0. 

Arguing exactly as in [7] Chapter 3] and [8] Chapters 3 and 4] , using also the argument in Appendix IDl 
and Footnote [2] we see that for every e £ (0, 1) there exists a connected open subset rig C £1 containing the 
chosen end of (M, g, k) and a solution u £ £ C^cC^o) °f tne J an g equation H(u £ ) + tr(k £ )(u e ) = on Hq with 
the following properties: 

(i) We have that {x £ : \x\ > A} C £Iq for some A > 1 that is independent of e £ (0, 1). The topological 
boundary 8Qq of fig is a smooth properly embedded hypersurface in whose components are either 
marginally inner trapped or marginally outer trapped with respect to the unit normal pointing into 
S7q. The components are A-minimizing with A := 1 + 2nsup ;Eg0 e e(o l) l^ £ ( a; )l m ^ ( m ^ ne language of 
[5]) and stable (in the sense of (fl9)) ). 

(ii) We have that u(x) — > as x — > oo in 57. We have that u e (x) diverges to plus infinity if x £ fig 
approaches a marginally inner trapped component of the boundary of S7g, and to minus infinity if 
x £ O converges to a marginally outer trapped boundary component. 

(iii) The graphs {(x,u £ (x)) : x £ 57g} of u £ are complete hypersurfaces of M x R that are stable and 
A-minimizing in f2. 

The A-minimizing property and stability of the graphs asserted in ([ml) implies that of the components of 
dflg in (P, cf. [30, p. 254], [7J Lemma A.l], and the discussion in Appendix [B] 

As in [5] or the discussion in Appcndix[B] we see that the stability and the almost minimizing property (via 
uniform local mass bounds) lead to curvature estimates for these graphs that are independent of e £ (0, 1). 
We now let e \ and pass the graphs of u £ to a smooth, properly embedded^ subsequential limit that 
contains a connected complete graphical component whose domain contains the asymptotically flat end and 
which satisfies all the above properties with e = 0. Using the mean value theorem and that there are no 
MOTSs or MITSs (with respect to k) in Cl besides 9Af_ and 8M + we see that this graphical component has 
all the properties asserted in the conclusion of the theorem. □ 

^The almost minimizing property does not imply uniform curvature estimates near dQ. As in Appcndix|B] we use the stability 
and the completeness for that. Once we know that a smooth limit exists, we can use the almost minimizing property to rule 
out sheeting. 
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Remark 3.2. Both in Proposition 13.11 above and in Theorem 13.31 below, the outermost property of M- and 
M_|_ prevents the domain of the sought-after graphical solution of the Jang equation from "popping outward" . 
If we think of the results in this section in analogy with the classical Jenkins-Serrin theory [16j H] > then this 
outermost property takes the place of the Jenkins Serrin flux conditions ((6j, (jJJ), and (jHJ in the statement 
of Theorem 16.11 below (with Hq = 0). 

Theorem 3.3. Let (M, g, k) be a complete asymptotically flat initial data set of dimension n where 3 < n < 7. 
Assume that dM + and dM- are both non-empty and that they intersect transversely. There exists a smooth 
solution u : M.nJlif — > K of the Jang equation H(u)+tr(k)(u) = such that u(x) — >■ as x oo in the end 
of (M,g, k) contained in M_nM + and such that u(x) —¥ ±oo as x € M—HM+ approaches an interior point y 
of a component of M^ndM± C d(M-C]M + ). The topological closure of the graph {(x,u(x)) : x £ M_nM + } 
of u in M x K is a smooth properly embedded hypersurface with manifold boundary (dM + fl dM-) x M. 

Proof. Let £1 := M_ n M+. We denote by v the unit normal vector field of the hypersurface M T n dM± 
pointing into f2. Let d-Q := M+ n dM- and d+Cl := M_ n 8M+. Note that cLfi and 9+f2 are manifolds 
with boundary and that their manifold boundaries coincide. The topological boundary of f2 is the disjoint 
union of <9_f2, d+fl, and dM- D dM + . Given a component 7 of d±£l, let 9 7 € C°°(7) be positive in the 
interior of 7 and zero on its manifold boundary. 

Let e £ (0, 1) be so small that the sets {exp y t<d 1 (y)i'(y) : y £ int 7 and t £ (0, 2e)} C f2 are disjoint as 
7 ranges over the components of d±Q. Let % e € C^ c (f2) with values in [—1, 1] be supported in the union of 
all these sets and such that x £ (y) = ±1 on Crf, := {exp y tQ 1 (y)v(y) : y £ int 7 and t € (0,e)} when 7 is a 
component of 9±fL 

Let 7 be a component of d-Vl. Consider the hypersurface {(exp y e(l— ■ e~ h )<d(y)is(y), h) : y £ int 7 and /i e 
(0,oo)} cOxi. It has piecewise smooth manifold boundary consisting of 7 x {0} and (dj) x [0,oo). It is 
the graph of a positive function tt? £ C^ c (Cr^). We have that H(u^) + tr(fc)(Sf,) = 0(e). 

Similarly, let 7 be a component of c? + 0. Consider the hypersurface {(exp y e(l — e )Q(y)v(y), h) : y £ 
int 7 and /i 6 (—00, 0)} C f2 x KL Its topological boundary is the union of 7 x {0} and (dj) x (—00, 0]. This 
hypersurface is the graph of a negative function u£ £ C^ c (Cr^). We have that H(yL) + tr(k)(yl) = 0(e). 

Choose c > so that for all e > sufficiently small H(vl) + tr(fc e )(il^) < — 2e on Cr^ when 7 is a 
component of <9_f2 and such that H{yV) + tr(k £ )(u^) > 2e on Cr^ when 7 is a component of <9+f2. Here, 
k E :— k + cex £ ■ 

Given t > 0, the functions — | + and | are, respectively, super and sub solutions of the regularized 
Jang equation H(u) + tr(fc e )(u) = t ■ u on Cr^. Let C > be a constant greater than HSup ee / n^gjj |fc e (x)|. 
Then y and — y are constant super and sub solutions of this equation. 

As in j7, Chapters 3 or 4] one sees that for every t > there exists a (Perron) solution uf £ C^ c (f2) of 
H(uf) + tr(fc e )(uj) = t ■ u\ such that — y < u\ < y on fl : such that u\ < — | + on Cr^ for components 
7 of c?_f7, such that | + < uf on Crf, for components 7 of and such that |itf(a?)| < &a(|^|) on 

{x £ fl : \x\ > A} for some A > 1 sufficiently large. Here, 6a is as in Appendix [l)R 

It follows that if y is an interior point of <9_f2, then limsupj.^ u\ {x) < — | and that if y is an interior 
point of then lim inf ^^.j, , x en u t ( x ) ^ f ■ 

The mean curvature of the graphs of u\ is bounded by 2C so that they are 2C-minimizing (in the language 
of [5]) in Q, x R, cf. [3 Appendix A]. A standard application of Allard's boundary regularity theorem exactly 
as in ,7, Chapter 4], using that the intersection dM- n 8M+ is transverse, shows that the closure of the 
graph of u\ in M x (-f , |) is a C 1,Q manifold with boundary (9M_ n dM + ) x (-f , f ). The C 1,Q estimates 
near the boundary depend only on C and the geometry of dM- n <9M + ; they are independent of e, t > 0. 



2 The results in [7] are for compact domains. One way to obtain u e t that decays to zero in the asymptotically flat end is as 
a limit of solutions u\' T on {x € Q : \x\ < r} with zero boundary values on {1 £ S] : |i| = r} as 1 < r -> 00. Note that 
|u^' r (a:)| < ft^d^l) on {x G n : A < |x| < r} for a fixed A > 1 sufficiently large. 
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We now pass the graphs of u\ to a geometric subsequential limit as e, t \ t 0. The existence and analysis 
of such limits is exactly as in [7J Chapters 3, 4] (which in turn are largely based on [30]). If in particular the 
limit along the subsequence (t n ,s n ) — > (0,0) were not a graph with the properties asserted in the theorem, 
there would be some x g CI such that the sequence {itj™ (x)}^ =1 is unbounded. For dehniteness, let us 
assume that u e t "{x) — > — oo, possibly after passing to a further subsequence. There exists a sequence of 
upward translations of the graphs of «f™ that converge, possibly after passing to a further subsequence, 
locally smoothly as hypersurfaces to a vertical cylinder Exl, where S C f2 is a smooth properly embedded 
submanifold with boundary cLM n d+M that encloses a bounded region CI with dCl with x £ Cl, and such 
that -He + trx;(fc) = 0. Here, the mean curvature is computed with respect to the unit normal pointing out 
of Cl. We can argue exactly as in (8[ Proposition 4.1 and Remark 4.1] that there exists a MOTS in M_ that 
is homologous to dM_ and which encloses {x} U M_. (The point is that we can force a blow down of the 
Jang equation in the complement of the closure of Cl in Af_.) This contradicts the assumption that 9Af_ is 
the outermost MOTS. □ 



4. Canonical blow up of the Jang equation 

In view of analogous results for Scherk-type minimal and constant mean curvature graphs on bounded 
domains, it is tempting to conjecture that the solutions of the Jang equation constructed in Proposition 
13.11 and Theorem 13.31 are unique with their properties. In Section [5] we prove such a uniqueness result in 
the special case where k = 0. In the case of general second fundamental form k, we will show in Theorem 
14.11 below that there exist canonical, pointwise maximal solutions of the Jang equation for example when 
M- C M + . The proof of this result proceeds via a curious geometric variant of the Perron method that uses 
the outermost condition built into the definition of M + in lieu of a super solution for the problem. The basic 
ingredients are variations of classical PDE techniques, cf. in particular [TBI [TjJ EU E2] and also [7J and the 
references therein, and the analysis of geometric limits of the Jang equation developed in [30, • 

Theorem 4.1. Let (M,g,k) be a complete asymptotically flat initial data set of dimension n, 3 < n < 7, 
fix an end, and let M + C M be the complement of the total inner trapped region of (M, g, k) with respect to 
that end. Assume that there exists a solution u : M + — » R of the Jang equation H(u) + tr(fc)(u) = such 
that u(x) — » oo as dist(ir, DM) — > and such that u(x) — > as x — > oo in the asymptotically flat end. The 
pointwise supremum of all such solutions is again a solution with the same properties. 

Proof. The results in Appendix ID] show that there exists A > 1 such that \u(x)\ < &a(|x|) on N := {x £ 
M + : \x\ > A} for every solution u of the Jang equation as in the statement of the theorem. 

Fix a smooth function u : M + — > M as in the statement of the theorem. Using translates of u to obtain a 
priori oscillation bounds and standard methods as in 7j Lemma 2.2], one sees that for every x G M + there 
exists p D (x) € (0, dist(a;, dM+)) small such that for every p £ (0,p D (x)) the equation H(«) + tr(k)(v) = 
on B p {x) with continuous boundary data on S p {x) admits a solution v £ C 00 (B p (x)) n C (B p (x)). The 
notion of (Perron) sub solutions u £ C(M+) and super solutions u £ C(M+) of the Jang equation H(u) + 
tr(k)(v) = on M+ can thus be defined in the usual way. Consider the class of functions S u = {u £ 
C(M+) : u is a Perron sub solution of the Jang equation, u > u on M+, and \u{x)\ < /3\(\x\) for all x £ 
M + with \x\ > A}. This class is closed under taking pointwise maximum and under lifting u £ S u to the 
functions £ C(M + ) that equals u on the complement of B p (x) and equals the solution v of H(v)+tr(k)(v) = 
on B p {x) such that !i = iion S p (x), for every r £ (0,p D (x)) and every x £ M + . Note that u £ S u . The 
function u p : M + ~ > RU{oo} is defined pointwise by u p {x) := sup ugSu u(x). Let fi := {x £ M + : u p {x) < 
oo}. Note that N C CI. 

We claim that CI is open, that u p \n is a smooth solution of the Jang equation, and that \im x ^,.y }X ^ u p (x) — 
oo for every y £ dCl. To see this, fix x £ CI. Following through the standard proof of the regularity of the 
Perron solution (cf. e.g. [13j p. 25]) we see that given p £ (0,p D (x)), there exist {u i }'^L l C <S„ such that 
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u < Ui < u.2 — ■ ■ ■ < u p on B p (x), such that H{u i ) + tr(/c)(u i ) = on B p (x) for all i = 1,2,..., and 
such that hnij-^oo ttj(x) = u p {x) < oo. The analysis of geometric limits of solutions of the Jang equation 
shows that the geometric limit of the graphs of in Be (x) x R contains the graph of a smooth solution 
u x ' r € C^ c (£Iq' p ) of the Jang equation above some open subset £l x ' p C B p (x). Moreover, lim^oo u^y) = oo 
for all y G Be[x) \ Qq' p , hm z _j. y , 2en ^.p u a: ' p (z) = oo for all y G S|(x) n <9^q' p , and Bp{x) n 9f2o''' i s a 
smooth properly embedded MITS in Be(x). (The point is that the functions are bounded below by u so 
that there can be no cylindrical components in their geometric limit, cf. the argument in Appendix lAl and 
the properties listed in Step 4 in Subsection 16.21 Note that f2g' p might have several components.) Clearly, 
u x ' p < u p on £1q' p . That u x,p — u p on the connected component of J1q' p containing x follows from the 
strong maximum principle for differences of solutions of the Jang equation, as in the standard proof of the 
regularity of Perron solutions. Since also u p > u everywhere on M+, we can deduce all the properties of 
and u p asserted at the beginning of this paragraph. 

The argument above shows that away from dM+ the boundary of f2 is a smooth properly embedded 
MITS. That the boundary of f2 is smooth and embedded up to d+M follows from the characterization in 
Appendix [B] of boundaries of domains that support solutions of prescribed mean curvature equations with 
infinite boundary data. The definition of M+ implies that fi = M+. Clearly, u p has all the properties 
asserted in the conclusion of the theorem. □ 

Remark 4.2. By taking the least super solution instead of the largest sub solution in the proof of Theorem 
14. II we obtain an analogous result where M+ is replaced by M_ and blow up to plus infinity at the boundary 
is replaced by blow down to minus infinity. 

5. Uniqueness of a blow up function u at the outermost minimal surface 

Let (M,g) be a complete asymptotically flat initial data set of dimension n, 2 < n < 7, with k = 0. Fix 
one of the ends and let £1 := M_ = M + be as in the introduction. In this case, d£l is called the horizon of 
(M,g). Note that dfl is a minimal surface. Let <9_f2 and be unions of different components of d£l such 
that dfl = d-flGd+n. 



The arguments proving Proposition ^ . 1 1 show that there exists a smooth solution u : ft — > R of the minimal 
surface equation 

\y/l +\DU\*J 

such that lim^-^y u(x) — ±oo for y £ d±Cl, and such that u(x) — > as |x| — > oo in Q. In this section we 
show that there is a unique function u with these properties. The proof is a straightforward adaption to our 
situation of a general argument due to J. Nitsche [26] as applied in e.g. [16j [17j [35j [15] to establish uniqueness 
of the Scherk-type graphs constructed there. We give the complete argument since the asymptotically flat 
ends require some care. 

To see that u is unique under the present assumptions, note first that the results in Appendix [B] show 
that the divergence of u near <9f2 is uniform in the distance to the respective components of the boundary, 
and that the upward and downward solutions of the graph converge geometrically to the vertical cylinders 
d + ft x R and cLfi x R respectively. In particular, 

Du 

(2) lim (x) = Tv(y) for y e d±tl, 

where v is the unit normal of dVl pointing into fl. Using the argument in Appendix [D] we see that 

(3) \u{x)\ + \x\\Du(x)\ = 0{\x\ 2 - p ) as |x| oo in fl 
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for every /3 € (2, n). 

Suppose that v : fi — > R is a second solution with these properties. Fix T € (0, oo) that is a regular value 
of both (u — v) and (v — u). Using ([2]) we see that 



sNOy{xefi:dist(x,an)=s,|(«-u)(x)|<T} \ \/l + \Du\ 2 + \Dv\ 

Here, rj denotes the unit normal of {x € SI : dist(x, <9f2) = s} pointing towards 90. Using the decay estimates 
(|3|) for u and v we obtain that 

lim [ {u- v )g(—£^ = -—£^ — ,^\dH n - x = Q 

r ^°° J {xen-.\x\=r,\u-v\<T} \^l + \Du\ 2 ^l + \Dv\ 2 J 

where r\ is the unit normal of {x € £1 : |x| = r} pointing towards the end. Using the divergence theorem and 
that u, v satisfy the minimal surface equation, this implies that 

0= / Jn„-n„—^ Dv Wn 



Using the strict convexity of the functions £ — ^ v/l + <?(x)(£, £) on T x Af for all x € M we conclude that the 
integrand is pointwise non-negative with equality at x € O if only if Du — Dv at x. It follows that u and v 
can only differ by a constant. Since we assume that they both tend to zero on the asymptotically flat end, 
we obtain that u = v, as desired. 

We see no way to extend this argument to non-zero A; at this point and have to contend with the existence 
of the canonical solution guaranteed under the hypotheses of Theorem 14.11 



Remark 5.1. Let u : fi — > R be as above. It follows that 

]x\)dn n -' 

Thus the unique blow up solution witnesses the area of the horizon at infinity. 

6. Extensions of the classical Jenkins-Serrin theory 



H n - l {d + n)-n n - 1 {d-n)= lim / g(Du,D\. 

r ^°° J{xm-.\x\=r} 



6.1. Introduction. The classical Jenkins-Serrin theory [TBJ [TT] characterizes those bounded domains £1 C 
R 2 with piecewise smooth boundary for which there exists a solution u : Q — > R of the minimal surface 
equation such that 

lim u(x, y) = ±oo and 

(x,y)-*-(xo,yo)Ed±Cl 

lim u(x,y) = <j)(x ,y ). 

(x,2/)->(xo,3/o)eoofi 

Here, the sets do^l, d-fl, and <9+f2 are unions of the smooth components of the boundary such that 
<9f2 = <9+£lU<9-ilU<9oOU{corners}, and cf> € C(doil) is a given function. A basic example of such a con- 
figuration is when ft = (— f , § ) x (— f , § ) and u(x,y) — log cosfaj • The g ra P n °f u m ^ 3 i s 01 course the 
classical Scherk surface. In general, these domains are precisely those for which the connected components 
of d±Vl are straight line segments such that no two segments in <9_f2 and no two segments in d + £l have 
an endpoint in common, for which the geodesic curvature of do^l is non-negative, and which satisfy the 
Jenkins-Serrin condition [T6l |T7] . When doil ^ these conditions demand that the circumference of every 
polygon inscribed in Q whose endpoints are chosen from the finitely many corner points is strictly greater 
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than twice the total length of its sides that coincide with segments in d+Cl and also greater than twice the 
total length of its sides that coincide with segments in d-tt. When do^l — the Jenkins-Serrin condition is 
the same except for the inscribed polygon that is the whole domain; one demands that the length of <9+fi 
equals the length of cLSl. 

An important and influential development of the field was accomplished by J. Spruck [35] , who has 
extended the classical Jenkins-Serrin theory to graphs u : O C R 2 -> R of constant mean curvature two. 
Provided that the piecewise smooth boundary dfl consists of a union <9 + f2 of circular arcs of unit radius 
that are convex towards f2, a union of circular arcs d-Q of unit radius that are concave towards f2, and a 
union of boundary arcs doQ whose geodesic curvature is greater than or equal to 1, he finds necessary and 
sufficient conditions for the existence of solutions u assuming (arbitrary) continuous boundary values on dafl 
and tending to oo on approach towards d+tt and — oo on approach towards d-fl provided a sub solution 



exists on the domain Q* obtained from flapping the negatively curved components of the boundary outward. 
(That this process gives a domain is a further additional assumption in [35].) The advantage of the piecewise 
smooth domain f2* over Q is that its boundary arcs are all convex, so that solutions of the constant mean 
curvature equation with prescribed boundary values away from the corners can be constructed on it. 

A further important contribution to the theory is due to U. Massari [20] . who has extended the Jenkins- 
Serrin theory to arbitrary dimensions and variable (Lipschitz) mean curvature in the case when d^tt =/= 0. 
His techniques arc different from J. Spruck's. In particular, no flapping of the boundary is required. In the 
special case where the mean curvature is constant, the necessary and sufficient conditions he provides for 
the existence of a solution are not obviously the same as those given in [55] . We note that condition (3.3) in 
[20] stands in, morally and technically, for the requirement ([4]) above that a sub solution of the prescribed 
mean curvature equation exists on the domain. In the book of E. Giusti |14| Chapter 16], an extension of 
the technique of U. Massari to the "slightly more complex" case where do ft = is presented in the minimal 
surface case. 

More recently, the Jenkins-Serrin theory for minimal graphs has been extended to H 2 x R in [25] and 
to M x R in [27j [22] (where (M, g) is a general complete Riemannian surface) . The Jenkins-Serrin-Spruck 
theory for constant mean curvature graphs has been developed for i 2 xl and S 2 x R in [T5J. A further 
important recent development are the results of P. Collin and H. Rosenberg [4] and A. Folha and S. Melo 
[10] who give necessary and sufficient conditions for the existence of Scherk-type minimal and constant mean 
curvature graphs on ideal polygons (with infinite area) in hyperbolic space. 

As an example of a particularly interesting application of the construction of Scherk-type graphs on Rie- 
mannian surfaces, we mention the surprising construction of harmonic diffeomorphisms between the complex 
plane and the hyperbolic space by P. Collin and H. Rosenberg [4]. 

In this section, we prove an extension of the Jenkins-Serrin-Spruck theory for domains CI C M with 
doCl — in Riemannian surfaces (M, g) and for general Hq £ [0, oo). When Hq = 0, our result is marginally 
different from the corresponding result in [57] because we consider the possibility of closed geodesies in the 
blow up/blow down analysis (in Step 4 below). In the case where Hq > 0, neither do we hypothesize the 
existence of a sub solution on fi, nor do we make any additional assumptions regarding the existence of an 



u : fl* -)• R 



(4) 
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auxiliary domain as in [351 115) . 

It will be apparent from our proof that our approach extends (for the most part) verbatim to appropriate 
domains in Riemannian manifolds of dimension 2 < n < 7. The description of admissible domains in higher 
dimensions is cumbersome. Also, the Jenkins-Serrin-Spruck conditions become impossible to verify in ex- 
amples that are not very symmetric, limiting the application of these results in higher dimension. For this 
reason, we omit a detailed discussion here. 

In our proof we work with solutions of the (finite boundary value) Dirichlet problem that will in general 
not assume any particular boundary data but only have the right asymptotic behavior. This is an important 
difference with the construction in [35j where solutions of the Dirichlet problem are constructed on an auxil- 
iary domains £1* whose boundary is sufficiently convex away from finitely many points to construct solutions 
that take on prescribed continuous data. The construction of fi* in [35] proceeds by flapping the negatively 
curved boundary components cLfl outward so they become convex arcs. This construction requires symme- 
try of the ambient manifold that is not in general available. For this reason, the construction in [35] [15] is 
carried out only in the simply connected space forms of dimension two. 



6.2. The case where dofl = 0. Let (M, g) be a complete boundaryless Riemannian surface, let Hq £ [0, oo), 
and let fi C M be a connected bounded open set such that <9fi = <9fi. Here and below, we will use "<9" to 
denote the topological boundary of a set. We assume that dfi is piecewise smooth and in fact the union 
of finitely many properly embedded arcs {Ai,Bj} and properly embedded closed curves {Ek,Fi} such that 
the outward geodesic curvature of each arc Ai and closed curve Ek is constant and equal to Hq and such 
that the outward geodesic curvature of each arc Bj and closed curve Fi is constant and equal to — Hq. We 
assume that all the curves and the interior of all the arcs are pairwise disjoint, and that no two arcs Ai and 
Aii and no two arcs Bj and By have an endpoint in common. 

The union of the closed curves {Ek} and the interiors of the positively curved arcs {Ai} is denoted by 
<9 + fi. The union of the closed curves {Fi} and the interiors of the negatively curved arcs {Bj} is denoted 
by cLfi. The endpoints of the arcs {Ai, Bj} are called the corners of <9fi. The assumptions imply (via the 
strong maximum principle when Hq > 0) that any two arcs that share an endpoint meet at a non-zero angle. 

A generalized polygon is a non-empty open subset P C fi with dP = dP and such that dP is piecewise 
smooth and consists of finitely many of the following building blocks: 

(1) Finitely many arcs of constant geodesic curvature Hq whose endpoints are amongst the corners of 
<9fi and whose interiors are embedded and pairwise disjoint. We also require that each arc whose 
interior intersects dfi is one of {Ai, Bj}. 

(2) Pairwise disjoint embedded closed curves of constant geodesic curvature ±Hq that either lie entirely 
in fi or coincide with one of {Ek,Fi} and which are disjoint from the boundary arcs. 

Theorem 6.1. Let (M,g) and fi C M be as above. A necessary and sufficient condition for the existence 
of a smooth function u : fi — > R such that 



(5) 
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is that 

(6) H 1 g (d + n)=H £l(n) + H 1 g (d_Q) 
and that 

(7) 2H 1 g {d+Q. n 8P) < U l g {dP) + H Q C 2 g {P) and 

(8) m x g (a_ n n ap) < w J (ap) - ff ^ (p) 

for every generalized polygon P C Q. 

The conditions ([6]) , (JTJ) , (|8]) appear in the classical work of Jenkins-Serrin (when Hq = and (M, g) is R 2 
with the Euclidean metric) and its generalization due to J. Spruck (when Hq > and (M, g) is Euclidean 
space), and also in the work of L. Hauswirth, H. Rosenberg, and J. Spruck [15] (where Ho > and (M,g) 
is one of R 2 ,S 2 ,IHI 2 with their constant curvature metrics). The necessity of these conditions follows from a 
standard argument, that we summarize briefly: 

Let u : 51 — >• R be as in the statement of Theorem 16.11 Let U C M be a non-empty and open subset 
such that U fl 5±51 = U n <951. The discussion in Appendix IBl shows that the graphs of the functions u =F t 
converge as hypersurfaces smoothly on compact subsets of U x R to <9±51 x 1 as f -} oo. In particular, 
as x G 51 approaches a point xo G c?±f2, the horizontal part of the downward unit normal of these graphs, 
X{x) := (1 + |_Dw| 2 ) _1 / 2 Z?w, converges to ± the outward pointing unit normal of c?±f2 at xq. The necessity 
of condition ([5]) follows from applying the divergence theorem to the vector field X on smooth interior ap- 
proximations of the domain 51. The necessity of conditions (J7J, © follows from the same argument applied 
to generalized polygons P C 51, using also that < 1 for x G dP n 57. 

Remark 6.2. A consequence of the the existence of a graph as in ([S]) is that every 7 G {Ai, Bj ,Ek,Fi} is stable 
in that J \Dtp\ 2 > J 7 (i? 2 + K g )(f) 2 for every i\> G C : (7) that vanishes on the boundary of 7 if 7 G {Ai, Bj}. 

Here, D is the (tangential) gradient and n g is half the scalar curvature of (M,g). This implies, for example, 
that in Euclidean space and when H Q = 1 , a domain that satisfies the Jenkins-Serrin-Spruck conditions ^ , 
(O, §8§ must also satisfy %\{Ai),'hL\{Bj) < it. The condition that %\{Bj) < it was part of the assumptions 
in [351 P- 16]. In fact, our proof (see property (jnj) in Step 4 below) shows that the conditions (O, (0), 
([5]) need only be verified for generalized polygons whose boundary arcs and closed curves are stable. This 
sharpening of the classical Jenkins-Serrin-Spruck condition is useful when constructing examples. 

Remark 6.3. The complement of a generalized polygon P in 51 is again a generalized polygon. Condition ((SJ) 
for a generalized polygon P C 51 follows from condition ([7]) applied to 51 \ P in view of ^ . 

Step 1: Construction of the auxiliary domain 51 

Fix a component 7 of <9-|-51 and a smooth function 7 G C°°{^) that is positive on 7 and which vanishes on 
its (manifold) boundary. Let v be the unit normal of 7 pointing out of 51. The piecewise smooth domain 51 
is obtained from 51 \ {vertices} by adding the crescents Cr 7 := {exp e (i9 7 (6')^(6')) : t G (0, e) and 9 G int(7)} 
as 7 ranges over all components of 5±51. Here, e > small is chosen so that there are no issues with the 
regularity of the exponential map and such that the crescents Cr 7 are pairwise disjoint. 

Step 2: Construction of barriers on Cr 7 and the functions H/.(x) 

Let 7 be a component of 5 + 51. We would like to find solutions that tend to 00 on approach to 7 and 
hence require a sub solution. The hypersurface {(exp e (ee h <d 1 (8)i'(9)), h) : h G (— 00, 0) and 6 G int(7)} of 
M x K is the graph of a (locally) smooth function u 7 : Cr 7 — > R whose downward unit normal corresponds 
to the outward unit normal of the cylinder. 
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Let 7 be a component of cLfi. We would like to find solutions that tend to — oo (and hence require a 
super solution). As above, the hypersurface {(exp e (ee~ h Q 1 (9)i'(9j),h) : h G (0, oo) and G intf/y)} is the 
vertical graph of a locally smooth function u 7 : Cr 7 — > M. 

if x € Q \ {corners}, 

if x G Cr 7 if 7 is a component of cLil, and 

if x G Cr 7 if 7 is a component of 9+51 

is locally Lipschitz. Moreover, H(x) = Hq + 0(e) uniformly on Cl. 

Let x : f2 — > [—1, 1] be a locally smooth function such that \ = il near Cr 7 for 7 G d±fl. Let k > 1. 
We define a locally Lipschitz function iJ^ on f2 by Hk{x) :— H{x) — fc _1 / 2 x(a;). Note that \fk + w 7 is a sub 
solution of the equation 

(9) div ( - DU 1 = H k + ~ ■ u 

on Cr 7 when 7 is a component of S+fi, and that — -v/fc + Sy is a super solution for this equation on Cr 7 when 
7 is a component of d-tt. 

Fix a constant C > sup fc>1 x£ f l \Hk(x)\. Then — Ck and Ck are, respectively, sub and super solutions for 
(J9j) on Cl. The introduction of the capillarity regularization (|9]) of the prescribed mean curvature equation 
so that large constants become barriers is exactly as in [30) . 

Step 3: The construction of 

Let Uk G C^™(i7) be the largest (Perron) sub solution of equation ([9]) that lies below Ck on all of Cl and 
below —\fk + u 7 on all crescents Cr 7 corresponding to components 7 of 9_f2. To justify the existence of 
such a solution, we refer to 32, p. 375], the interior gradient estimate stated in Appendix |A1 and also |36l 
Theorems 1.1 and 1.4], [131 Chapter 16], and [3 Chapter 3]. The maximum principle implies that Uk lies 
above — Ck on all of Cl and above ^fk + u^ on all crescents Cr 7 corresponding to components 7 of d+Q. From 
(j9|) we see that the mean curvature of the graph of Uk is bounded uniformly by 2C on Q. 



The function H : — >• R defined as 
H(x) 




Step 4: Geometric limits of graph(ufc) 

We claim that there is a subsequence of {uk} and there exist disjoint open subsets CIq, f2- of Cl 

and u G C it j"(fio) with the following properties: 

(a) Cl = (Hq Uf!_ U 0+) n Cl. In particular, the topological boundaries <9f2o, c?Sl_ , 951+ of Hq, f2 + locally 
separate fio, 57 + from their respective complements f2\f2o, 57\J7 + in Cl. Moreover, dfloHCl, 
dfl- fl f2, and c?Sl + n Cl are properly embedded C 2,a hypersurfaces in Cl. 

(b) For every x G f2+ there exists an open neighborhood of x in f2 + so that (y) exceeds a given constant 
for all y in this neighborhood, provided i is sufficiently large. Put differently, diverges to plus infinity 
locally uniformly on Q + . 

(c) For every x G 0_ there exists an open neighborhood of x in f2_ so that (y) lies below a given constant 
for all points y in this neighborhood, provided i is sufficiently large. Put differently, Uk { diverges to minus 
infinity locally uniformly on fi_. 



14 



MICHAEL EICHMAIR & JAN METZGER 



(d) We have that — > u in Cf^(flo). In particular, 

,. / Du 



= H on 57q. 



\y/l + \Du\- 
(c) Cr 7 C Cl± when 7 is a component of <9±57. 

(f) The sets 957on(57U {corners}), S57_ n (57 U {corners}), and <957_ (1 (57U {corners}) consist of finitely many 
arcs and closed curves of constant geodesic curvature in 57. These arcs and closed curves are pairwise 
disjoint in 57. The arcs are properly immersed and embedded in 57. Their endpoints are corners of 57, and 
the endpoints of any one arc may coincide. The closed curves are contained in 51 and they are properly 
embedded. 

(g) If 7 is a component of <957-|- n 357o, then lim xe n ,:E-Kc u(x) — ±00 uniformly near xq € int(7). 

(h) The geodesic curvature of a component 7 of <957 H 951+ is constant and equal to H when we orient 7 by 
the unit normal v pointing into 57 + . Every divergent series of downward translations of the hypersurface 
graph(u) = {(x,u(x)) : x <E 57o} converges to (dQo H 957+) x t in C 2 ' a on compact subsets of 57 x R. 

(i) The geodesic curvature of a component 7 of 957o fl <957 _ is constant and equal to Ho when we orient 7 
by the unit normal v pointing into Slo- Every divergent series of upward translations of the hypersurface 
graph(u) = {(x,u(x)) : x g 57o} converges to (dflo H 957_) x 1 in C 2,a on compact subsets of 57 x R. 

(j) The geodesic curvature of a component 7 of (957_ n c>57_|_ is constant and equal to Hq when we orient 7 
by the unit normal v pointing into 57+. There exists a cylindrical neighborhood of 7 x R in 57 x R in 
which the graphs {(x,Uki(x)) : x € 57} converge in C 2,a to 7 x R on compact subsets. 

(k) The graphs {(ar, Uk^x)) : x € 57} converge as embedded C 2,a hypersurfaces on compact subsets of 57 x R 
to the union of the cylinders (57 n c?o57) x R, (9_57 n (9+57) x R and the graph {(x, u(x)) : x € 57o}, as 
i — > 00. 

(1) The vector fields ^7^=^ | 2 are locally equicontinuous in 57. 
(m) With 7 and v as in (|h| - |j]) we have that 

In fact, the integrand on the left converges to 1 locally uniformly on int(7). 
(n) The arcs and closed curves 7 in (jQ, (0), and those in Q interior to 57 are stable in the sense that 

f (H 2 + K g )i) 2 <M} g < f \D^\ 2 dHl for all ip G C 1 ( 7 ) with supp(0) C int( 7 ). 

Here, k 9 is half the scalar curvature of (M, g) and D is the (tangential) gradient of ip. 
The properties listed here extend classical results about limits of monotone sequences in the Jenkins-Serrin- 
Spruck theory, cf. [16j [T71 [35l [15l [27]. For limits of general (i.e. not necessarily monotone) sequences, 
some of these properties can be inferred directly from the results of L. Mazet [23]. The ideas in [53] have 
been employed to prove Jenkins-Serrin-Spruck type results for minimal graphs supported on domains in 
Riemannian surfaces in [2 2) and for constant mean curvature graphs in [10j . 

The above properties are also variations of classical results on generalized solutions of the minimal surface 
equation |21) or the geometric theory of Jang equation [30] . Below, we supply a few comments on the proofs 
of these properties to assist the reader. The proofs extend to domains of dimension n < 7. 

The particular barriers used in the construction of Uk in Step 4 imply (gj). 



Properties (jaj) - (jmj) can be deduced from the compactness and regularity properties of almost minimiz- 
ing boundaries, of which graphs of bounded mean curvature are a basic example, along with the geometric 



JENKINS-SERRIN TYPE RESULTS FOR THE JANG EQUATION 



15 



Harnack principle (as in Appendix [XJ , which ensures that geometric limits of our graphs are made up of 
graphical and cylindrical components. For limits of minimal graphs, this is the approach of |21) . For the case 
of geometric limits of solutions of the regularized Jang equation (including the capillarity regularization) , 
this approach has been worked out in detail in [7] , to which we refer the reader for more general statements 
statements and further references. 

Property (JTJ) follows from Lemma IC.ll in Appendix [Cj 

The argument leading to (jnj) is very similar to that in Appendix [Bj The Jacobi identity (fl"7f is replaced 
by the differential inequality 

.3 , nu. |2 , Tj„ f, ,,„:< . . ..I n, it. 1 ... , I ■ ,. > -> n, . 



A Gfc ^ + (l^fcl 2 + Rc 5+(dx 3)2(n,^))^ - -9 ^ 1 =±== ,D(H k + ~ • u k )\ 4 < k-^lDxWl on Q 

where all geometric quantities on the left are computed for the graph Gk := {(x,Uk(x)) : x G Q}. The addi- 
tional contribution to the stability inequality (|18p disappears when we take geometric subsequential limits of 
Gk and its vertical translates as n — > oo. Since the arcs 7 for which (jnj is asserted appear as cross-sections 
of vertical cylinders that appear in such limits, and because \h\ 2 + Rcg+^dx 3 ) 2 {v, v) reduces to Hq + n g on 
such cross-sections, we are done. 

Step 4: Analysis of the limit using the Jenkins Serrin Spruck conditions 

The analysis of the geometric limit of the graphs of the solutions : il — > K using the Jenkins-Serrin- 
Spruck condition here shares structural features with the analysis in [35j (see in particular Sections 5 and 6 
therein) or [151 Section 7], Because of the absence of a sub solution in particular, our main technical step, 
Case b below, is quite different. 

The properties listed in Step 4 show that the components P of fio n O, fi + fl 0, and f2_ D O are generalized 
polygons in Q. If P is a component of fi_ fl f2, then 

(10) H C 2 g (P) + Urn sup / ^-dC] > 

i— >oo Jp "*i 

H 1 g (dPnd + n)+'H 1 g (dPnn) + iwiiv£ [ j^'f"^ dn i 

Here, v is the unit normal pointing out of ft. Note that the second term on the left is always non-positive. 
Similarly, if P is a component of f2+ fl then 

(11) H C 2 g (P)+1imM [ ^dC 2 g < 

% — y 00 J p YZ^ 

-Hl(dPnd-n)-'H 1 g (dPnn) + ]imsu.p [ 9^,Du ki ) , 

1^00 JdPnd+n yl -I- \Duk z \ 2 
where again the unit normal v points out of O. The second term on the left is always non-negative. 

We analyze the following cases: 

Case a: ^ n Q C Q 

Let P be a component of f2_ fl f2. The Jenkins-Serrin-Spruck condition for P implies that 

h £ 2 (p) < Hl(dPnd + n) + Hl(dPnsi)-HUdPnd-Si) 
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which leads to an immediate contradiction with (|10[) . since clearly 

<ttj(apn0_fi). 



9i ^ Du ^ -dHl 



dPnd-n \J\ + \Duki\ 



Thus Case a cannot occur. 



Case a': jt 0+ n ft C ft 

Let P be a component of il + n SI. The Jenkins-Serrin-Spruck condition for P implies that 

h c 2 g {p) > n\{dP n d+n) - «*(0P n ft) - ^(cip n d_ft) 

which leads to an immediate contradiction with since clearly 

g{y, Du ki 



dH 1 

dPnd+n V 1 + \Duki\ 2 



< Hl(dPnd + n). 



Thus Case a' cannot occur. 



Case b: SJcfi- 

The Jenkins-Serrin-Spruck condition for P = ft_ n ft = ft implies that 

JJ o £2(fi)=w£(0 + fi)-Wj(0_fi). 

In conjunction with (|10[) we conclude that 

limsup [ ^dC 2 a = and that liminf [ 9("> Du i°i) m i = _ u lr g _™ 

Passing to a further subsequence, if necessary, we see that 

(12) lim £ 2 ({x e ft : ^ < -£» = for every e > 
and, using (JTJ) , that 

(13) lim ' = — 1 locally uniformly on <9_ft 

^/l + \Du ki \ 2 

where v is the unit normal pointing out of ft. 

Fix a component 7 of cLft and let z e 7. It follows from the assumptions that z G ft_. Consider the 
functions u ki {x) :— Uk i (x) — u ki (z). (This is an upward translation for fcj large.) Then 

A . ( _Duki \ rr , . u kz {z) 
div — = = /ft. + — . 

\y/l + \Du ki \*J ki ki 

Using that ^^(a;)] < Cki for all x € ft we see that the mean curvature of these graphs is uniformly bounded. 
We pass to a further subsequence so that k~ l u ki (z) converges to a constant c G [— C, 0]. We pass to a further 
subsequence so that the graphs of the u ki converge geometrically in C 2 a to a union of properly embedded 
graphs and cylinders on compact subsets in ft x R. The mean curvature of these graphs and cylinders at 
a point (x,x 3 ) € ft x R in the geometric limit is H{x) + c. The point (z, 0) € 7 x R is contained in the 
geometric limit. Using ([T^| we see that (z, 0) is contained in a cylindrical component 7 x R of the limit, 
where 7 C O is a properly embedded curve whose mean curvature at a; € 7 is given by H(x) + c. Moreover, 
the tangent spaces of 7 and 7 agree together with their orientation at any point of 7 n 7. 
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We claim that 7 = 7. In particular, c = 0. To see this, we distinguish two cases. 

First, assume that 7 C Cr 7 . Then the assertion is a consequence of the maximum principle. (We use that 
c < here.) 

Second, assume that 7 n SI 7^ 0. Recall that every y £ 7 has an open neighborhood in SI that is separated 
by 7 into two components such that tends to plus infinity locally uniformly in one component as k — > 00, 
and such that tends to minus infinity locally uniformly in the other component. Since 7 fl SI ^ 0, we 
conclude that {x £ SI : limsup^^ fc~ Uk ( (x) < 0} = {x £ SI : lira sup^^ k^ Uk i (x) < c} contains a 
non-empty open subset. In conjunction with (fT2|) we conclude that c = 0. It follows that in this case, 7 
and 7 satisfy the same geometric equation. Further, we know that they intersect non-trivially, and that at 
any point of intersection they intersect tangentially with the same orientation. The Hopf boundary point 
lemma shows that 7 = 7, which contradicts the assumption that 71"! SI ^ 0. Hence 7 = 7 is the only possibility. 

The argument in the preceding paragraph shows that there exists a relatively open neighborhood U 7 of 
7 in SI that is disjoint from the crescents corresponding to the positively curved boundary components such 
that converges to — 00 locally uniformly in f7 7 fl Cr 7 , and to 00 in J7 7 fl SI. 

We can repeat the above reasoning for any of the components 71,..., 7 m of cLSl, choosing a point 
zi £ ji for each component. Passing to a further subsequence and relabeling if necessary, we may assume 
that Uki(z) > Uk^Zi) for all i £ {l,...,m} where z = z\. Also, we may pick y £ SI near 71 so that 
> Uki f } y ^ > > Uh i — > 0. Finally, let := Ukf — Ufe 4 (j/). For k large, this is an upward translation. Then 

divf - °^ 

For the remainder of the argument, we replace Cr 7 by t/ 7 (1 Cr 7 (keeping the same notation) for all com- 
ponents 7 of d-Sl. This may shrink the auxiliary domain SI slightly. However, we gain that Uki(%) ~> — 00 
locally uniformly in these new crescents Cr 7 . Note that — > 00 locally uniformly in Cr 7 when 7 is a 
component of 9+S1. We can take a subsequential geometric limit of the graphs of 11^ just as we did for the 
original sequence it^ so that (jaj) - (JTJ) continue to hold. We use SIq, Sl_, Sl + instead of Slo, SI-, Sl + to avoid 
confusion. As before, the components of Sin Slo and SlnSl± are generalized polygons. We claim that SI C Slo- 
To see this, note that (y, 0) is contained in the geometric limit. This implies that y £ <9S1_ U dSl + U Slo so 
that SI c Sl± is impossible. The cases ^ SI n Ct± C SI can be ruled out exactly as in Cases a and a' above. 
It follows that SI = Slo, and we can conclude as in Case c below. 

Case b': SI C Sl + 

Exactly as in Case b, we can conclude that a sequence of downward translations of Uk t will converge to 
a solution of the original problem. We point out that in fact, the analysis can be shortened considerably in 
this case because large constants are super solutions of the equation. 

Case c: Slo = SI 

In this case, the solutions it^ converge to the sought-after solution u : SI — > K by jsj) and (jgj). 

Remark 6.4. Let (M, g) be a complete Riemannian manifold of dimension n with 2 < n < 7. Let H £ C°°(M) 
and let SI C M be a bounded domain whose boundary can be written as the disjoint union of hypersurfaces 
cLSl and 9+ SI such that Hq_^{x) > H(x) with respect to the unit normal pointing into SI and such that 
Hd + n{x) < H(x) with respect to the unit normal pointing out of SI. There exists an open subset U C SI 



j =H k . I " fc - l Ufcl( 
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containing a neighborhood of d-fl whose boundary in Q is a smooth hypersurface £ whose mean curvature 
at x £ £ with respect to the unit normal pointing out of U equals H(x). Such a set U can be found by 
minimizing the brane functional 

(14) U ^H"-\nnd*U)- f H{x)dC n {x). 

This was proven by M. Fuchs in [IT] Theorems 2.1 and 4.1]. The existence of a hypersurface with prescribed 
mean curvature H also follows from the non- variational approach in [3] [7] ; simply note that such surfaces 
are MOTSs in the initial data set (M,g, k := — ^zjg). The proofs in [3j[7] proceed by constructing a limit of 
solutions of regularized Jang equations whose boundary values diverge to plus and minus infinity near <9+f2 
and <9_f2 respectively. The observation we have exploited in the proof of Theorem 16.11 that the capillarity 
term in the regularized Jang equation contributes "with a good sign" in the flux integrals (fT0|) and (1TTT) can 
be used to show that the boundaries of prescribed mean curvature arising in this way also minimize (|14p . 
Cf. Remark 3.2]. 

Appendix A. Remark on the Interior gradient estimate for the prescribed mean curvature 

EQUATION 

Let (M,g) be a complete Riemannian manifold, let O C M be a non-empty open subset, and let u £ 

loc 



C,'"(Cl) be a solution of the prescribed mean curvature equation 



(15) div f ■ ° - 1 = H + 1 ■ u on Q 

\^l + \Du?) 

where t € (0, 1) and where H is a bounded, locally Lipschitz function. The geometric operators appearing in 
this equation are with respect to the metric g, so that the left hand side represents the scalar mean curvature 
of the graph G(u) := {(x,u(x)) : x € 0} c Q x M with respect to the product metric g + (dx n+1 ) 2 and 
the downward pointing unit normal. Then, given x G fl, the gradient of u at x is bounded in terms of the 
geometry of (M, g) near x, the size of H in Lipschitz norm near x, and a bound on u near x. This interior 
gradient estimate is derived in |33j using iteration techniques, and e.g. in |18j . |19j . [3 6) using the maximum 
principle in a clever way. The estimates for the gradient obtained in these references are explicit. 

Here, we include a short, indirect, and conceptually simple proof of this well-known interior gradient esti- 
mate. This proof works when the dimension n of M lies in the range 2 < n < 7. It is based on the regularity 
theory of (almost) minimizing boundaries. The ingredients of the proof are all classical. For convenience, 
our references below are to the recent paper [7] whose notation and language are compatible with ours here. 
There, the reader will find detailed references to the literature. 

Suppose the estimate fails (with a counterexample in dimension 2 < n < 7). Then there exist (M,g), 
57 C M, igS), and H as in the statement, a precompact open subset B C f2 containing x, and a sequence of 
solutions Ufc € C, of (fT5"j) with t = tk£ (0, 1) such that Uk{x) — 0, such that |ttfc(j/)| < T for some T > 
and all y £ B, and such that \Duk{x)\ — > oo as i — > oo. Let Gk denote the graph of Uk above B, and let 
= (1+ l-Dwfcl 2 ) -1 / 2 denote the vertical component of its upward pointing unit normal. By assumption, 
the mean curvatures of the graphs Gk are uniformly bounded independently of i. Hence the graphs Gk are all 
contained in a precompact class of C 1,a almost minimizing (relative) boundaries inBxR (cf. [7j Appendix 
A and Remark 4.1]). Moreover, there exists a constant (3 > depending only on the Lipschitz norm of H 
on B such that A Gfc ^ +1 < /3i/™ +1 (cf. p3 (5) and (6)]) holds weakly on G k . Moreover, (x,0) € G k for 
every i. Let G be a subsequential varifold limit of G k as k — > oo. Then G is a properly embedded C 2,a 
two-sided hypersurface in B x R, and the subsequence of the G k approaches G in C 2,a on compact subsets 
of B x R. Let v n+1 denote the vertical component of the "upward" unit normal of G (the orientation is 
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inherited from Gk)- Then v n+1 is non-negative and Agi/" +1 < f3v n+1 weakly on G. By the strong maximum 
principle, on any component of G, v n+1 is either everywhere positive (so that the component is a graph) or 
everywhere vanishing (so that the component is a vertical cylinder). See [2TJ Section 4.2] and Lemma 
2.3] for references and a more detailed discussion of this geometric Harnack principle. Since \Duk{x)\ — > oo, 
it follows that the connected component of G containing (x. 0) is cylindrical. This is impossible, since all 
the graphs Gk and consequently also G are contained in the slab M x [— T, T]. 

Appendix B. Domains of solutions of prescribed mean curvature equations with infinite 

BOUNDARY DATA 

Let (M, g) be a complete Riemannian manifold of dimension n, 2 < n < 7, let ^ 57 C M be a non-empty 
proper open subset, and let u € C (57) be a solution of the equation 



Du 



y/T+\Dv\- 



where Hq G R is a constant. Let U C M be an open set such that U fl 957 ^ and assume that 
lim^-^y^gn u(x) = oo for every y £ U (1 951, where 957 is the topological boundary of 57 in M. 

The topological boundary U fl 957 of 57 relative to U has the structure of a smooth two-sided properly 
immersed constant mean curvature hypersurface i : £ — > U. More precisely, given y 6 J7 fl 957, there exists 
an open neighborhood O y of y in U and a diffeomorphism ^ : O v — > i?™~ 1 (0) x (—1, 1) with 0y(y) = such 
that one of the following conditions hold: 

(1) There exists a smooth function / : B™ _:L (0) -> K with /(0) = 0, d/(0) = 0, and \df{x)\ < 1 for all 
x e B^- l (0) such that <j> y (ClnO v ) = {(x,t) e B^-\0) x (-1,1) : t < f {x)}. 

(2) There exist smooth functions f : £^ -1 (Q) R with f{0) = 0, df l (0) = 0, |df*(a;)j < 1, and 
fHx) < P{x) for alia; E B™ _1 (0) and i € {1,2} such that <j> y {nnO y ) = {(x,t) e Bf _1 (0) X (-1, 1) : 
* </ 1 M or/ 2 (x) <t}. 

Moreover, for every tp € C^(S) we have that 

(16) / (N 2 +Rc 9 (^))V 2 < / \D^\l 

Here, h is the second fundamental form and v is the unit normal vector field of the immersion i : S — > J7, 
integration is with respect to the volume form of the pull-back metric g = l* g, and D is the gradient with 
respect to g. 

For n = 2, the conclusions here are (roughly) [16l p. 329] (for Hq — 0) and [32l Theorems 6.1 and 6.2] (for 
Hq ^ 0). For n > 2, it is shown in Chapter 8 [35] that U D 957 is a hypersurface of constant mean curvature 
Hq assuming that 957 is a C hypersurface. 

The proof of the general result presented here may be well-known to experts. For completeness, we include 
a sketch of how these statements may be derived conveniently from techniques applied in [3 [8] . 

Note that the graph of u, G := {(x, u(x)) : x € U D 57}, is a properly embedded hypersurface in U x R. 
Let V <e U be a non-empty open set with smooth boundary that is compactly contained in U. Using the 
divergence theorem for the vector field (1 + \Du\ 2 )~ 1 / 2 (—Du, 1) in the region {(x, x n+1 ) : x G 57 and u(x) < 
x n+1 } n (V x [5, T]) for regular values S, T of u with S < T, we obtain that 

H 9 V +1)2 (Gn(y x [S,T])) < 2C n g {V) + |H |(T- S^-^dV). 
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In particular, we obtain locally uniform area bounds for G and its vertical translates. 

Let v n+1 := (1 + |Z3u| 2 ) -1 / 2 denote the vertical component of the (upward) unit normal vector field 
v = (1 + |Z)u| 2 ) _1 / 2 (— Du, 1) of G. The Jacobi identity (expressing the fact that vertical translation of G 
leaves the mean curvature constant) implies that 

(17) A G (v n+1 ) + (\h\ 2 + Rc g+(dxn+lv (v,is)y i+1 =0. 

Here, Aq is the (non-positive) Laplace-Beltrami operator of G with respect to its induced metric g. Let 
tp € C\{U x M) be a test function, multiply ([17]) by ^> 2 (z/ n+1 ) _1 , integrate over G, integrate by parts, and 
use the elementary estimate — g(D(t() 2 (i' n+1 )^ 1 ), Dv n+1 ) > — \Dip 2 \g to obtain that 

(18) f (\h\ 2 +Rc g+{dxn+1) 2(v, V ))i; 2 < f \D^\l 

JG JG 

The stability-based regularity theory of [29] applies and provides curvature estimates for G. In fact, the 
curvature of G n (V x [T, T + 1]) is bounded for every V <s [7 independently of T e 1. (The case n = 7 
requires an additional argument, cf. Remark 4.1 in [7].) Let Gt := {(i,utW) : x € 0} with ut '■= u — T. 
Note that these hypersurfaces are naturally ordered and that their area and curvature is bounded indepen- 
dently of T on compact subsets of U xl. Hence we can pass them to a geometric (varifold) limit as T /*• oo. 
It is elementary to check that the support of this limit equals (U n dfl) x M. The arguments in Appendix 
A of [SJ, in particular Corollary A.l and Remark A. 3, explain carefully how the asserted structure of the 
cross-section of this geometric limit follows from this, in particular why no more than two sheets can come 
together. That the stability property (fT8|) passes to the cross-section in the form (|16l) is checked exactly as 
in [301 P- 254]. 

Similar arguments characterize the boundary of open subsets fl of a Riemannian manifold (M, g) that 
support solutions u : — > R with infinite boundary values on U (1 dQ of the prescribed mean curvature 
equation 

div ( DU \ - F ( — \ 

™ \y/l + \Du\*J ~ [ X, y/1 + \Du\z) " 

Here, F is a smooth function on M x S n ~ 1 (M), where S n_1 (M) is the unit sphere bundle of (M,g). The 
boundary of Q in U then has the structure of a two-sided immersion i : S — > U such hat Hs(i(a)) = 
F(i(a),v(a)), where v is the unit normal field along the immersion pointing out of fi. Instead of (|16l) one 
obtains that for every e > there exists a constant C £ = (e, | Rcm |, |i?|c0 such that 

(19) (1-e) / |/i|V 2 < / \Dip\l + C e [ for every V e Cj(S). 
See (3)] for details. 

Appendix C. Equicontinuous vector fields from solutions of the prescribed mean 

CURVATURE EQUATION 

The following lemma is an immediate consequence of the classical compactness and regularity properties 
of graphs of bounded mean curvature in low dimensions, see e.g. [7j Appendix A and Remark 4.1] for precise 
statements and references. Though particularly useful in the context of Jenkins-Serrin type problems, it 
does not seem to appear explicitly in the literature. 
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Lemma C.l. Let (M,g) be a complete n- dimensional Riemannian manifold, 2 < n < 7 , let Q C M be a 
non-empty open subset, and let C > 0. The collection of continuously differ entiable vector fields 



J- := < — = : u £ C 2 (Q) and 



,. / Du 

div 



< c 



is equicontinuous on compact subsets o/O. 

Appendix D. Barriers for the Jang equation near infinity 

In (301 p. 248], certain rotationally symmetric barriers for the Jang equation were constructed on large 
subsets of the asymptotically flat ends of three dimensional initial data sets. The following proposition, 
which we quote from [6], is a straightforward extension of the construction in |30j to higher dimensions. 

Proposition D.l. Fix (3 £ (2, n). For A > 1 define 

f°° ds 
B\(r) := A / — on [A,oo). 

Then /3a is continuous and positive, it is smooth on (A, oo), and we have that ^ (r) tends to — oo as r \, A. 
There exists a constant c = c(/3) > 1 such that b\(r) < cA(r/A) 2- ^ and such that 6a (A) > c _1 A. 

Let (M,g, k) be an asymptotically flat initial data set of dimension n, n > 3, such that g^kij = 0(|x| - ' 3 ) 
as x — > oo in the asymptotically flat ends^There exists Aq — Aq(M, g, k, (3) > 1 such that for every A > Aq 
we have that H(b A (\x\)) + tr(fe)(6 A (|a:|)) > and H(-b A (\x\)) + tr(k)(-b A (\x\)) < on {x £ M : \x\ > A}. 

In this paper, we repeatedly use the barriers from Proposition ID. II along with the maximum principle, 
exactly as in [501 p. 249], to conclude that a solution u of the regularized Jang equation H(u) + tr(fc)(u) = 
t ■ u (with t > 0) that is defined on the complement of a compact subset of (M,g, k) and which tends to 
zero at infinity actually lies between — 6 A (|a;|) and 6a(|x|) provided that A > Ao is large enough so that 
{x £ M : \x\ > A} is contained in the domain of u. 
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